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Abstract. We introduce a real-space renormalisation group procedure for 
driven diffusive systems which predicts both steady state and dynamic properties. 
We apply the method to the boundary driven asymmetric simple exclusion process 
and recover exact results for the steady state phase diagram, as well as the 
crossovers in the relaxation dynamics for each phase. 



1. Introduction 

Nonequilibrium phenomena are observed throughout nature and society, and can often 
be thought of as some kind of driven diffusive system ^ . Examples include chemical 
kinetics, the motion of interacting molecular motors along a substrate, population 
dynamics in an ecology, the exchange of wealth or commodities in an economy, traffic 
flow and many more. The ongoing interest in driven systems has been sustained by 
the wide variety of behaviour the model systems exhibit, which incorporate not only 
elements familiar from equilibrium statistical mechanics, such as phase transitions, 
cooperative behaviour, and so on, but also many dramatic new features. These include 
phase transitions induced by external driving and even in one dimension; long-range 
(i.e. power-law) correlations, not only at criticality, but generically in the steady state; 
and very rich dynamics. 

Model driven diffusive systems are typically defined on a lattice on which particles 
hop from site to site, where the precise definition of the stochastic particle dynamics 
is motivated on physical grounds. Nonequilibrium steady states are constructed by 
driving a current of particles through the system. Despite being simply stated, they 
still exhibit the wide range of phenomena. Indeed, a virtue of their simplicity is that 
in some cases the models can be solved exactly. However, most cases elude exact 
solution; moreover, even in cases where the steady state can be written down, the 
task of computing macroscopic observables may be too difficult. For these reasons, it 
is necessary to develop approximate techniques. 

A technique that has been hugely successful in equilibrium statistical mechanics 
is the renormalisation group procedure. Here, we show how one can adapt 
these equilibrium methods to nonequilibrium systems, where we devise a way to 
accommodate the feature that distinguishes the nonequilibrium steady state from its 
equilibrium counterpart — the existence of currents — in the renormalisation group 
prescription. In particular, our approach includes the effects of fluctuations, which, 
in certain equilibrium and nonequilibrium systems, destroy the transitions predicted 
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and removal from the right boundary 



by mean field theory. Thus our approach builds on previous prescriptions (2j, and it 
complements real-space techniques devised for quantum-spin chains, and adapted to 
stochastic models |3|E]- 

We test our approach using the asymmetric simple exclusion process with open 
boundaries, which we introduce and review in section |21 In section |21 we outline 
the matrix product formulation of the master equation. This formulation provides 
a convenient notation for the purposes of the steady state renormalisation which we 
describe in section^] We recover the exact critical point and find qualitative agreement 
for the phase diagram. We extend the renormalisation to the dynamics in section|31and 
again we obtain qualitative agreement with the exact predictions for the relaxation 
dynamics. 

2. The asymmetric simple exclusion process (ASEP) 

The asymmetric simple exclusion process with open boundaries is one of the most 
widely studied driven diffusive systems. Exact results show that it undergoes both 
first- and second order phase transitions in the steady state and exhibits different 
relaxation dynamics in each of its phases. For these reasons, it represents an ideal test 
bed with which to assess the efficacy of approximate methods. 

The model is defined on a chain of L sites, labelled I = 1, . . . , L. The occupation 
variable n; = if site / is vacant and n/ = 1 if it is occupied; multiple occupancy is 
forbidden. The stochastic dynamics are illustrated in figure ^ in continuous time, 
a particle at site I < L hops to its nearest neighbour site to the right with rate p, 
provided it is empty. If site 1 is vacant, it becomes occupied with rate a, and if there 
is a particle at site L, it leaves the system with rate /3. 

The phase diagram for this model is known exactly 0. In the limit of large L, 
depending on the values of the ratios a/p and f3/p, the system may reside in one 
of three possible phases (see figure for a/p < 1/2 and a < (3, the system is in 
a low density phase with particle density p < 1/2; for f3/p < 1/2 and a > (3, the 
system is in a high density phase, p > 1/2; for a/p > 1/2 and (3/p > 1/2 the current, 
J, through the system assumes its maximum value J = 1/4 and the density p = 1/2. 
The transition between the high and low density phases, indicated by the dotted line 
in the phase diagram, is first order — along the transition line a low density region 
coexists with a high density region and the system density p — 1/2. The transition 
between these low current phases and the maximum current phase, indicated by the 
solid line in the phase diagram, is second order. 

The relaxation dynamics have also been obtained recently using the Bethe ansatz 
: in the maximum current phase the relaxation time scales as where the dynamic 
exponent z = 3/2; along the coexistence line z = 2; in the two low current phases 
away from the coexistence line, the relaxation times are finite. 
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Figure 2. Exact phase diagram for the boundary driven ASEP. 



3. Matrix product formulation of the ASEP 

The exact solution for the steady state of the ASEP was obtained using a matrix 
product technique 0- In this approach, the aim is to represent the probability that 
one observes the system in a configuration n = (rti, . . . , n^) at time t, P{n, t), as 

L 

P{n, t) = Z-\W\ ^{niD + (1 - ni)E)\V) , (1) 
i=\ 

where D and E are (in general time-dependent) matrices, and and \y) are (time- 
independent) vectors. Here, Z is a normalisation given by 

Z = {W\C^\V). (2) 

where we have introduced the (time-independent) matrix C = D + E. Thus, a string 
of matrices is used to represent a configuration of the particles in the system, where 
the matrix D is used to represent an occupied site and E a vacant site. 

3.1. Steady state algebra 

In the limit i ^ oo, it was shown in that and (O do indeed give the steady 

state probabilities provided the matrices D and E and the vectors {W\ and \ V) satisfy 
the algebra 

D + E = DE , (3) 

{W\E = [p/a){W\ , (4) 

D\V) = {p/l3)\V) . (5) 

Moreover, matrix representations were found which satisfy (O-lO- We remark however 
that in order to compute the probabilities to observe a particular configuration, it is 
not necessary to find such matrix representations; it is instead sufficient to use ^ 
to order all the matrices E to the left of the matrices D, and then use the boundary 
conditions and © to obtain probabilities in terms of a, /3 and p. 

3.2. Dynamic algebra 

In [7] , it was shown how the steady state algebra can be generalised to a full dynamic 
algebra such that (Q) and (0) represent the probabilities at all times t. In this case. 
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Figure 3. Blocking the lattice into groups of b = 3 sites. 



the matrices and vectors must satisfy 



D=[C-\K] , 



(6) 
(7) 
(8) 
(9) 



AC-^D = DC-^A , 
{W\{a/p~C-'^A) =0 
{(3/p-AC-'\V)^0, 



where A = DE is a current operator (capturing the fact that in order for a particle to 
hop there must be a particle to the left of a vacancy). 

4. Steady state renormalisation of the ASEP 

The key feature which distinguishes nonequilibrium steady states from equilibrium 
ones is the presence of currents. Moreover, in driven diffusive systems with open 
boundaries, the boundary dynamics determine the particle current and density 
throughout the bulk of the system. This is in marked contrast to equilibrium systems 
where boundary effects are typically localised to the vicinity of the boundary. 

In this section, we present a renormalisation group method which accounts for 
these currents by employing a blocking procedure designed to accommodate the 
influence of the boundaries throughout the bulk of the system. 

4.1. Blocking transformation 

A blocking transformation is implemented by grouping sites of the lattice into blocks 
where each block contains b sites, as illustrated for the case 6 = 3 in figure |21 These 
blocks are then identified with the sites of a renormalised lattice with a number of 
sites L' = L/b. This must be supplemented by a prescription which identifies particle 
configurations within a block with a single coarse-grained block variable. We achieve 
this by using a majority rule identification, which for 6 = 3 is 



The issue for any such renormalisation group treatment then is what quantities 
are held fixed under the blocking and which are rescaled. We keep the current across 
each bond J and the bulk density p fixed, and rescale the parameters a/p and (3 /p. 
To achieve this, we consider pairs of adjacent blocks separately, and replace the effect 
of particle exchange with the rest of the system with particle injection and ejection at 
the boundaries of the pair of blocks, with rates a/p and (i/p, as illustrated in figure 
0] Thus the current and density in the rescaled system are written 



D' = DDD + DDE + DED + EDD , 
E' = EEE + EED + EDE + DEE . 



(10) 
(11) 
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Figure 4. Blocking transformation at fixed current for the boundary driven 
ASEP with b = 3. 



{W'\D'E'\V') 

{W'\C'^\V') ' ^ ' 

^ {WWW) 

^ {w'\C'\v') ■ ^ ' 

Under the majority rule matching between configurations, pU|) . the block current 
J and the block density p assume the form 

(W\{D^+D^E+DED+ED^){E^+E^D+EDE+DE^)\V) 

Jw\cW) 

_ {W\{D^+D^E+DED+ED^)\V) 

W\cW) ■ ^ ^ 

Note that we impose that neither the vectors and {W\ and \V) nor the matrix algebra 
Q to (O are changed under the rescaling. This reflects the fact that both the 
scaled and the unsealed configurations are taken to evolve under the same stochastic 
dynamics: any extension of the parameter space is prohibited under the rescaling. 
Finally, the rescaling of the parameters a/p and fi/p is achieved by setting J' = J 
and p' — p. 



4.2. a = 13 

Up to this point, the use of the matrix product formulation has been nothing more 
than a choice of notation. The task now is to evaluate the correlation functions (|12|l 
to IjlSfl and this is relatively straightforward if one makes use of the algebra |j3Jl to 
((SJ, though we remark that even here the algebra only simplifies the calculations: 
one could evaluate the relevant correlations directly from the master equation, for 
example. For the case a = j3 one obtains (still for & = 3) 

p' = p = 1/2 , (16) 
which is also the result obtained in the exact solution for a = (3, and 
2x' _ 66a; + 83a;2 + 46a;3 + 12x4 

3a;'2 + 2x' ~ Mx + 126a;2 + 112a;3 + TOx* + aOa;^ + ' ^ ' 

from matching J' = J, where x = p/ot. Stable fixed points of this equation, 
{a/p)* = {a/p) = {a/p)', are found at {a/p)* = and {a/p)* = 3.959, and these 
are separated by an unstable fixed point at {a/p)* = 1/2. At the unstable fixed point, 
J* = 1/4. Thus the {a/p)* = fixed point represents the low-current phase, and the 
transition between this phase and the maximum current phase is given by the unstable 
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3.248 
0.423 


0.423 
3.248 


3.959 
3.959 


0.415 

oo 


oo 
0.415 


2.795 

oo 


oo 
2.795 



Table 1. Fixed points of the steady state blocking transformation for 6 = 3. The 
values given as decimals are evaluated numerically — all the rest are determined 
analytically. The fixed points depicted in the flow diagram shown in figure |5] are 
labelled A to G. 
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Figure 5. Schematic flow diagram in the a/p-f3/p plane for generic values of b. 
The locations of the the fixed points A to D do not depend on b, the remaining 
fixed point values are fe-dependent. 

fixed point, for which we recover the exact value (a/p)* — 1/2. Moreover the current 
and density at this point are also given by their exactly known values throughout the 
maximum current phase: J* = 1/4 and p* = 1/2. 

4.3. a^P 

For a ^ P the scaling equations are straightforward to obtain, but rather lengthy so 
we do not include them here. Their fixed point structure is given in table ^ It is 
also straightforward to obtain the fiow diagram by making some choice for an initial 
set of parameters and seeing how they evolve under several iterations of the scaling 
transformation. The result is given in figure^ and should be compared with the exact 
phase diagram shown in figure [21 

The fully unstable fixed point A controls the second-order transition between the 
low current and maximum current phases. At this point, the current and density 
assume their exact critical values {pc, Jc) = {\, \)- The flow line AE is a separatrix 
which marks the boundary between the high density regime in the low current phase 
(represented by point B) and the maximum current phase. Similarly, the line AF 
separates the low density regime in the low current phase (represented by point C) 
from the maximum current phase. Points B, C and D are all zero current fixed points. 
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The first-order transition in the model is represented by the flow line AD, which 
separates (exactly) the high density from the low density regimes in the low current 
phase. Also, across the line a/p + (3/p = 1, the gradient of the density profile near 
the boundary changes sign. This is the line BAG in the flow diagram. 

Thus we have interpretations for the fixed points A to D. For the remaining fixed 
points, we have no such physical interpretation, though their presence is necessary 
for the flow diagram to adopt the same qualitative structure as the phase diagram. 
Critical exponents can be obtained by linearising the scaling equations in the vicinity 
of the critical fixed point. Also, the scaling outlined above was illustrated for blocking 
six sites to two {h — 3). It is also possible to perform the rescaling for eight sites to 
two [b — 4) or nine sites to three (5 = 3), for example, thereby increasing the block 
size or the size of the renormalised clusters. Though these considerations do not effect 
the topology of the phase diagram, they may improve predictions for exponents. We 
illustrate this point more fully, for dynamic scaling, in the next section. 



5. Dynamic scaling for the ASEP 

So far, the scaling has concerned only the statics. In this section, we present a blocking 
transformation in order to learn about the relaxation to the steady state. The blocking 
proceeds in the same spirit as in the previous section and exploits the dynamic algebra 
© to @. Then the scaling transformation, linearised in the vicinity of the physical 
steady state fixed points (A to D in the fiow diagram) , informs us of the nature of the 
relaxation dynamics in each of the corresponding phases. 



5.1. Blocking transformation 

We begin by replacing blocks containing two sites with a single renormalised site (and 
will generahse subsequently). This single site is parameterised by {a/p)' and [P/p)' 
and its dynamics are governed by the time evolution of {D'{t)): 

djP'jt)) ^ {W'\D'it)\V') 

dt {w'\C'\v') ' ^ ' 

(c.f. equation H13|l ). Again, the coarse-grained block variable is identified with 
configurations in the original system through a majority rule statement which, for 
blocks containing two sites, assumes the form 

D' = DD + \{DE + ED)^\{DC + CD) . (19) 

This, carried into l|18|l . yields the evolution in the blocked system: 

d{D'{t)) _ (^1 [D{t)C + cm] \V) ^^^^ 



dt 2{W\C'^\V) 

Now the dynamic algebra is used to compute the time evolution of the relevant 
correlation functions. Although the rhs of (|20|l depends only on {DC) and {CD), 
the evolution of these two correlation functions depends on the third time-dependent 
function {DD). Thus the time evolution of the two-site block can be expressed as a 
matrix equation, where the relaxation time is determined by the largest eigenvalue of 
a 3 X 3 matrix. This eigenvalue gives the gap to the steady state eigenvalue (which is 
zero), and we refer to this gap as A*^^'=), where Nc represents the number of sites in 
the block: here, Nr — 2. 
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Now, under the coarse-graining, we hold the gap fixed therefore we set 

AiK) ^ A(We) ^ (21) 

where A*^^':-' is the gap in the blocked system and N'^ is the number of sites in the 
coarse-grained system: here, N'^ — 1. Thus the dilation factor b = N^/Nc- Through 
equation 1)21(1 we obtain a scaling equation in terms of the boundary parameters. 
Because only a single quantity — the gap — is fixed under scaling, we can derive a 
scaling relationship for a single parameter only. Thus we consider the case Up = (3p. 
This is sufficient to probe the dynamics near the steady state fixed points A and D. 
The scaling relationship between {a/p)' and a/p is not in general linear, therefore it 
is possible that the dynamics at the zero-current fixed point and those at the critical 
fixed point are characterised by different values of the dynamic exponent z. Linearising 
the scaling equation near either of these fixed points obtains 

= Aap, (22) 

which for each fixed point specifies the associated eigenvalue and exponent A = 6^^. 



5.2. a ^(3 

We now illustrate this procedure for Nc = 2 and iV^ = 1 in more detail, in the case 
a — f3. The dynamic algebra to © applied to equation lfTH|) yields 

{D{t))^{a/pY~2{a/pY{Dit)}, (23) 

and the gap is easily seen to be A'^^ = 2{a/p)' where the prime is used to indicate 
a parameter in the coarse-grained system. For the two-site system, the dynamic 
algebra applied to (|20fl leads to the evolution of the correlation functions which can 
be represented as 

/ {DD) \ ( -2a/p «/p \ / {DD) \ 

dt {DC) U 1 + {DC) ,(24) 

V {CD) I \ -I 1 -a/p I V {CD) I 

from which the gap is found to be A^^^ = a/p. Thus we obtain the scaling equation 

= Yp ' 

therefore A = 1/2 so z = 1. We remark again that one could just as easily perform 
this scaling working directly with the master equation: the operator algebra simply 
provides a neat formulation for our purposes. 

In this case the same dynamic exponent is obtained for zero-current and critical 
fixed points, but this is peculiar to blocking from two sites to one. In other cases, one 
obtains for {a/p)' a complicated nonlinear function of a/p. This function is linearised 
in the vicinity of the steady state fixed points A and D in order to compute the dynamic 
exponent in the phases they represent. The results of such a linearisation, for several 
values of Nc and A'^^ are shown in table El We also show the exact values which have 
only recently been obtained [H]. Our results appear to capture the crossover from 
z = 3/2 dynamics in the maximum current phase (represented by point A) to z = 2 
dynamics along the coexistence line (represented by point D). 
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z at a/p = P/p = 


z at a/p = 13 Ip =1/2 


2 ^ 1 


1 


1 


3 ^ 1 


1.118 


1.072 


3^2 


1.319 


1.196 


4 ^ 1 


1.194 


1.111 


4^2 


1.388 


1.221 


4^3 


1.486 


1.257 


5 ^ 1 


1.249 


1.133 


5-^2 


1.437 


1.234 


5^3 


1.531 


1.265 


5^4 


1.589 


1.274 


Exact 


2 


3/2 



Table 2. Estimates of the dynamic exponent near the zero current and J = 1/4 
fixed points (points A and D in Figure IS! . 



5.3. 13^1- a 

One can also obtain scaling equations for the case /3 = 1 — a in order to probe the 
dynamics at points B and C in the flow diagram, points which represent the high and 
low density phases, respectively. One finds that at both B and C, {a/p)' = {ct/p) 
for all Nc and N'^. This implies there is no system size dependence in the leading 
contribution to the gap for large L, which is consistent with the exact prediction of 
finite relaxation times in the low current phases away from the coexistence line 

6. Conclusion 

We have developed a renormalisation group procedure which we have applied to 
the boundary driven asymmetric simple exclusion process. In the steady state 
transformation, exact results for the critical point are recovered as well as a flow 
diagram that shares the qualitative structure of the exact phase diagram. The dynamic 
scaling also captures the crossovers in the exactly computed relaxation dynamics in 
each phase of the ASEP. 

These scaling methods seem to be widely applicable: they can be generalised to 
higher dimensions, multiple particle species and multiple site occupancy. Since only 
small clusters give, through the scaling, even the highly correlated properties, the 
scaling does not depend on the availability of exact results, and may be particularly 
useful in cases where the exact steady state can be written down but any calculation 
of macroscopic properties remains difficult. They are also direct and straightforward 
to interpret, without requiring any translation from quantum spins or abstract fields 
as in field-theoretic methods. 
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